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Abstract -It was recently found that cascading failures can cause the abrupt breakdown of a 
system of interdependent networks. Using the percolation method developed for single clustered 
networks by Newman [Phys. Rev. Lett. 103, 058701 (2009)], we develop an analytical method for 
studying how clustering within the networks of a system of interdependent networks affects the 
system's robustness. We find that clustering significantly increases the vulnerability of the system, 
which is represented by the increased value of the percolation threshold p c in interdependent 
networks. 
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Introduction. — In a system of interdependent net- 
works, the functioning of nodes in one network is depen- 
dent upon the functioning of nodes in other networks of 
. the system. The failure of nodes in one network can cause 
nodes in other networks to fail, which in turn can cause 
further damage to the first network, leading to cascad- 
ing failures and catastrophic consequences. Power black- 
outs across entire countries have been caused by cascad- 
ing failures between the interdependent communication 
and power grid systems [Ud]. Because infrastructures in 
our modern society are becoming increasingly interdepen- 
dent, understanding how systemic robustness is affected 
by these interdependences is essential if we are to design 
infrastructures that are resilient [3H6]. In addition to re- 
search carried out on specific systems [7HT6] , a mathemat- 
ical framework [17] and its generalizations [T8H20] have 
been developed recently. These studies use a percolation 
approach to analyze a system of two or more interdepen- 
dent networks subject to cascading failure [2TJ[22] . It was 
found that interdependent networks arc significantly more 
vulnerable than their stand-alone counterparts. The dy- 
namics of cascading failure are strongly affected by the 
structure patterns of network components and by the in- 
teraction between networks. This research has focused 
almost exclusively on random interdependent networks in 



which clustering within component networks is small or 
approaches zero. Clustering quantifies the propensity for 
two neighbors of the same vertex to also be neighbors of 
each other, forming triangle-shaped configurations in the 
network [23H25] . Unlike random networks in which there 
is very little or no clustering, real- world networks exhibit 
significant clustering. Recent studies have shown that, for 
single networks, both bond percolation and site percola- 
tion in clustered networks have higher epidemic thresholds 
compared to the unclustered networks [26H3Tj . 

Here we present a mathematical framework for under- 
standing how the robustness of interdependent networks is 
affected by clustering within the network components. We 
extend the percolation method developed by Newman [26] 
for single clustered networks to coupled clustered net- 
works. We find that interdependent networks that exhibit 
significant clustering are more vulnerable to random node 
failure than networks without significant clustering. We 
are able to simplify our interdependent networks model — 
without losing its general applicability — by reducing its 
size to two networks, A and B, each having the same num- 
ber of nodes N. The N nodes in A and B have bidirec- 
tional dependency links to each other, establishing a one- 
to-one correspondence. Thus the functioning of a node in 
network A depends on the functioning of the correspond- 
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ing node in network B and vice versa. Each network is with d! = s' — n' ■ The generating function of P' s , t , is 
defined by a joint distribution P st (generating function 

G {x,y) = J2Tt=o Pstx s y t ) that specifies the fraction of G (x,y,p) 
nodes connected to s single edges and t triangles [25] . The 
conventional degree of each node is thus k = s + 2t. The 
clustering coefficient c is 
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Therefore, the generating function of the remaining net- 
work after attack is 



Site Percolation of single clustered networks. — 

We begin by studying the generating function of remain- 
ing nodes after a fraction of (1 — p) nodes is randomly 
removed from one clustered network. After the nodes 
are removed, we define t\ to be the number of triangles 
of which node i is a part, d[ to be the number of single 
edges that form triangles prior to attack, and n£ to be the 
number of stand-alone single edges prior to attack. This 
network is thus defined by the joint distribution P n ',t',d'- 
The probability that a node has n' single edges from sin- 
gle edges is the sum of all the probabilities that nodes 
with more than n' single edges will have exactly n' edges 

°° ( s \ _ ' 

remaining, which is Q\(n') = X) / ) P n (1 ~ pY ™ ■ 

s=n' \ 12 J 

Similarly, the probability that a node has t' triangles is 
the sum of all the probabilities that nodes with more than 
t' triangles will have exactly t' triangles remaining. Since 
the probability that a triangle will survive is p 2 , the sum 
t 



is Q 2 (f) = £ 
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The probability 



that a triangle corner will have one edge broken is 1 _ 
and the probability that it will have both edges broken 
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Thus the probability that a node had d! sin- 
gle edges forming triangles prior to their destruction is 

t-f \ r 2p(l-p) irf' r (1-p) 2 U-t'-d' 



Combining 



these three, we have the corresponding generating function 
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G (xp + 1 - p, yp 2 + 2zp{\ - p) + (1 - pf). (2) 



We define s' = n' + d! to be the total number of single 
links of a node after attack. The joint degree distribution 



after attack is P' s , t , which satisfies P' s , v — Y^,n'=o Pn', 



G (x, y,p) =G (xp + l- p, yp 2 + 2xp(l -p) + (1 - p) 2 ). 

(4) 

The size of the giant component g(p) of the remaining 
network according to Ref. [26] is 



where 



g(p) = 1 - G (u,w 2 ,p), 

U = G q (u,V 2 ,p), 

v = G r (u,v 2 ,p), 
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and G q (x,y,p) = ± 9Go< £ y ' p) , G r (x,y,p) _ 
where and v are the average number of single links and 
triangles per node, respectively. 

As an example, consider the case when (1— p) fraction of 
nodes are removed randomly from a network with doubly 
Poisson degree distribution 



s! t!' 
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where the parameters fi and v are the average numbers of 
single edges and triangles per vertex, respectively. Accord- 
ing to Eq. ([TJ, the clustering coefficient is c = 2 v+(p+2v)' i • 
Then, G (x,y) = e K^-i) e Hy-i) an d G {x,y,p) = 
G q (x,y,p) = G r (x,y,p) = e lw+Mi-p)»](*-i) e »p 2 (y-i) ^ 
and u = v = 1 — g(p), leading to 

g(p) = 1 - e 



[fj.p+2p(l-p)v]g{p) e vp (g(p) 



-2s(p)) 
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This equation is a closed-form solution for the giant com- 
ponent g(p) and can be solved numerically. The criti- 
cal case appears when the derivatives of the both sides 
of Eq. (JSJ arc equal. That leads to the critical condition 
(k}p c = 1, which is independent of clustering. However the 
degree distribution of the doubly Poisson model changes as 
we keep the average degree and change the clustering coef- 
ficient. When the degree distribution is fixed, the critical 
threshold actually increases as clustering increases [29f30] . 
Furthermore, Fig. Q] shows the resulting giant component 
as a function of p. Note that single networks with higher 
clustering have smaller giant components. 
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Degree-Degree Correlation. — When constructing 
clustering in a network, it is usually impossible to avoid 
generating degree-degree correlations. To better under- 
stand the effect of clustering on degree-degree correlations, 
we present an analytical expression of degree correlation 
as a function of the clustering coefficient for a doubly 
Poisson-clustered network — see Eq. ([7]) . 

The degree-degree correlation [32] can be expressed as 
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where N m is the total number of m hop walks between all 
possible node pairs including cases i = j. 

The generating function of the degree of a node in 

oo 

the network is E Pstz s+2t = Gq(z,z 2 ). Let q st be 

s,t=0 

the fraction of nodes with s single edges and t trian- 
gles that are reached by traversing a random single link, 
where s includes the traversed link and r st is the frac- 
tion of nodes with s single edges and t triangles reached 

sP 

by traversing a link of a triangle, q s t = Vt , r s t 
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Their corresponding generating functions are G q (x,y) 
^^ilx and G r (x, y) = Moreover, N 3 = 

J2J2 a ij^U)' wn( 3re Nz(j) is the total number of two- 

i j 

hop walks starting from node j. The number of three- 
hop walks from a node i is equal to the total number of 
two-hop walks starting from all of its neighbors. Thus, 
N 3 = ^2 kjN2(j), where the number of two-hop walks 

j 

starting from a node j with degree kj will be counted kj 
times in JV 3 . Equivalently, N 3 = NJ2(s + 2t)P s . t N 2 (s,t), 

St 

where A^(s, i) is the number of two hop walks from 
a node with s single edges and t triangles. The gen- 
erating function of the number of single edges and of 
triangles reached in two hops from a random node is 
G 2 (x,y) = £P M • G q (x,y) ■ Gf{x,y). The generat- 

St 

ing function of the total number of links and of trian- 
gles reached within three hops starting from all nodes is 
G 3 (x,y) = N^Ps.t ■ (G q (x,y)) s{s+2t) ■ {G r (x,y)) 2t(s+2t) . 

St 

The number of fc-hop walks can be approximated by 
its mean in a large network 
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When both s and t follow a Poisson distribution, 
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which together with Eq. © leads to 
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where c is the clustering coefficient, Eq. (jTJ) . 

Figure [2] shows the relation between the degree correla- 
tion and the clustering coefficient c for a Poissonian net- 
work [see Eq. (0], for two given average degrees ((fc) =3 
and 4). The figure shows a positive degree-degree corre- 
lation across the entire range, which means the model is 
assortative [29]. The degree-degree correlation increases 
until c achieves half of its maximum and then decreases 
to zero when c reaches its maximum. When c is or the 
maximum, the nodes connect to either all single links or 
all triangles, respectively. 

Percolation on Interdependent Clustered Net- 
works. — To study how clustering within interdepen- 
dent networks affects a system's robustness, we apply the 
interdependent networks framework [17j . In interdepen- 
dent networks A and B, a fraction (1 — p) of nodes is first 
removed from network A. Then the size of the giant com- 
ponents of networks A and B in each cascading failure 
step is defined to be p\, p n , which arc calculated 

iterativcly 



p n = |U„_i5^(/i„_i),n is odd, 
Pn = p„gB(p n ), n is even, 



(11) 



where /iq 
isfy 



p and fi n are intermediate variables that sat- 



Pn =P3A(A*n-i),n is odd, 
Pn = pgs(A*n-i),n is even. 



(12) 



As interdependent networks A and B form a stable 
mutually-connected giant component, n — > oo and fx n = 
/i ra -2, the fraction of nodes left in the giant component is 
Poo. This system satisfies 



x =P9a(v), 
V =P9b(x), 



(13) 



where the two unknown variables x and y can be used 
to calculate p^ = xgs(x) = ygA{y)- Eliminating y from 
these equations, we obtain a single equation 



x = pgA[pgB{x)]. 



(14) 
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The critical case (p = p c ) emerges when both sides of this 
equation have equal derivatives, 

1 = P 2 ^^[p9b{x)]^-(x)\ x=Xc ^ Pc , (15) 

which, together with Eq. ([T4"]). yields the solution for p c 
and the critical size of the giant mutually-connected com- 
ponent, Poo(p c ) = x c g B {x c ). 

Consider for example the case in which each network has 
doubly-Poisson degree distributions as in Eq. (JT)). From 
Eq. (|13p . we have x = p{\ — Ua), y = p(l — Ub), where 

ua = va = e^ Ay+2y ^ 1 ^ v ^ A ^ UA ^ +VA ' p2 ^ A ~ 1 \ 

If the two networks have the same clustering, /i = ha = 
Hb and v = Va = vb, Poo is then 

Poo =p{l-e v ^-^ +2v ^f. (16) 

The giant component, Poo, for interdependent clustered 
networks can thus be obtained by solving Eq. ([TBI) . Note 
that when v = we obtain from Eq. (Til)]) the result ob- 
tained in Ref. [17] for random interdependent ER net- 
works. Figure I3a[ using numerical simulation, compares 
the size of the giant component after n stages of cascading 
failure with the theoretical prediction of Eq. (fTT|) . When 
p = 0.7 and p = 0.64, which are not near the critical 
threshold (p c = 0.6609), the agreement with simulation is 
perfect. Below and near the critical threshold, the simu- 
lation initially agrees with the theoretical prediction but 
then deviates for large n due to the random fluctuations of 
structure in different realizations [T7] . By solving Eq. (fT5|) , 
we have p^ as a function of p in Fig. [3b] for a given aver- 
age degree and several values of clustering coefficients and 
in Fig. [4a] for a given clustering and for different average 
degree values. As the figure shows, when higher clustering 
within a network is introduced, the percolation transition 
yields a higher value of p c (see inset of Fig. 13b]) . 

When clustering changes in this doubly Poisson distri- 
bution model, degree distribution and degree-degree cor- 
relation also change. First, to address the influence of 
the degree distribution, we study the critical thresholds of 
shuffled clustered networks. Shuffled clustered networks 
have neither clustering nor degree-degree distribution but 
keep the same degree distribution as the original clustered 
networks. The brown dashed curve in Fig. I3bl represents 
the giant component of interdependent shuffled clustered 
networks with original clustering c = 0.2. The figure shows 
that the difference in p c between the c = network and the 
shuffled c = 0.2 network is only 0.01, while the difference 
between the c = and the c = 0.2 networks is 0.12. In 
addtion, c = 0.2 clustered networks has no degree-degree 
correlation (Fig. [2]), which means the 0.12 shift of p c is 
due to clustering and not to a change in degree distribu- 
tion. We also show the critical thresholds of interdepen- 
dent shuffled clustered networks as the red dashed line in 



the inset of Fig. I3b[ Note that the change of degree distri- 
bution barely shifts the critical threshold. We next discuss 
the effect of the degree-degree correlation on the change of 
critical threshold. From Ref. [34] , the degree assortativity 
alone monotonously [??? JDM] increases the percolation 
critical threshold of interdependent networks. Because in 
our case degree-degree correlation first increases and then 
decreases (see Fig. [2]), while critical the threshold of in- 
terdependent networks increases monotonously [???] as 
clustering increases, we conclude that clustering alone in- 
creases the value of p c . Thus clustering within networks 
reduces the robustness of interdependent networks. This 
probably occurs because clustered networks contain some 
links in triangles that do not contribute to the giant com- 
ponent, and in each stage of cascading failure the giant 
component will be smaller than in the unclustered case. 

We also study the effect of the mean degree (k) on the 
percolation critical point. Figures Ha] and [4b] both show 
that, when clustering is fixed, the percolation critical point 
of interdependent networks decreases as the average degree 
(k) of network increases, making the system more robust. 
Figure I4bl also shows that a larger minimum average de- 
gree is needed to maintain the network against collapse 
without any node removal as clustering increases. 

Conclusion and Summary. To conclude, based 
on Newman's single network clustering model, we present 
a generating-function formalism solution for site perco- 
lation on both single and interdependent clustered net- 
works. We also derive an analytical expression, Eq. ([TO]) , 
for degree-degree correlation as a function of the cluster- 
ing coefficient for a doubly-Poisson network. Our results 
help us better understand the effect of clustering on the 
percolation of interdependent networks. We discuss the 
influence of a change of degree distribution and the degree- 
degree correlation associated with clustering in the model 
on the critical threshold of interdependent networks and 
conclude that p c for interdependent networks increases 
when networks are more highly clustered. 
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Fig. 1: Size of giant component g(p) in single networks 
with degree distribution Eq. (JT)) and average degree (k) = 
4, as a function of p, the fraction of remaining nodes after 
random removal of nodes. Curves are from theory Eq. [8l 
symbols are from simulation. 
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Fig. 2: Degree-degrees correlation as a function of the 
clustering coefficient for Poisson network (Eq. ([7])) with 
average degree (fc) = 3 and 4. Curves are from theory 
(Eq. I10[) and symbols from simulations. 
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Fig. 3: (a) Size of mutually connected giant component 
as a function of cascading failure steps n. Results are for 
c = 1, p = 0.64 (below p c ), p = 0.66 (at p c ) and p = 0.7 
(above p c ). Lines represent theory (Eqs. (fTTj) and (O) 
and dots are from simulations. Note that at p c there are 
large fluctuations, (b) Size of giant component, p^,, in 
interdependent networks with both networks having clus- 
tering via degree distribution Eq. ((7|) and average degree 
(fc) = 4, as a function of p. Dashed lines are number 
of interactions (NOI) before cascading failure stops ob- 
tained by simulation [33]. The star curve is for shuffled 
c = 0.2 network, which keeps the same degree distribution 
but without clustering and without degree-degree correla- 
tion. Inset: Green squares and solid line represents critical 
thresholds, p Cl of interdependent networks as a function of 
clustering coefficient c. Red dashed line represents critical 
threshold of shuffled interdependent networks which orig- 
inally has clustering coefficient c. The shuffled networks 
have zero clustering and degree-degree correlation, but has 
the same degree distribution as the original clustered net- 
works. In all figures, symbols and dashed lines represent 
simulation, solid curves represent theoretical results. 
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Fig. 4: (a) Size of giant component as a function of p for 
fixed clustering coefficient c = 0.1 and different average 
degrees. From right to left (k) = 3,4,5, ...,9. (b) Critical 
threshold p c as a function of average degree for different 
clustering coefficients. The solid curves are for interdepen- 
dent networks and the dashed curve is for single networks. 
Symbols and curves represent simulation and theoretical 
predictions respectively. 
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